Search for: All records

Abstract Haystack and Owens Valley Radio Observatory observations recently revealed strong, intermittent, sinusoidal total flux-density variations that maintained their coherence between 1975 and 2021 in the blazar PKS 2131−021 (z= 1.283). This was interpreted as possible evidence of a supermassive black hole binary (SMBHB). Extended observations through 2023 show a coherence over 47.9 yr, with an observed periodP_{15 GHz}= (1739.8 ± 17.4) days. We reject, withp-value = 2.09 × 10⁻⁷, the hypothesis that the variations are due to random fluctuations in the red noise tail of the power spectral density. There is clearly a physical phenomenon in PKS 2131−021 producing coherent sinusoidal flux-density variations. We find the coherent sinusoidal intensity variations extend from below 2.7 GHz to optical frequencies, from which we derive an observed periodP_optical= (1764 ± 36) days. Across this broad frequency range, there is a smoothly varying monotonic phase shift in the sinusoidal variations with frequency. Hints of periodic variations are also observed atγ-ray energies. The importance of well-vetted SMBHB candidates to searches for gravitational waves is pointed out. We estimate the fraction of blazars that are SMBHB candidates to be >1 in 100. Thus, monitoring programs covering tens of thousands of blazars could discover hundreds of SMBHB candidates. 

Abstract This paper addresses, for the first time, a key aspect of the phenomenology of compact symmetric objects (CSOs): the characteristics of their radio spectra. We present a radio-spectrum description of a complete sample of high-luminosity CSOs (CSO-2s), which shows that they exhibit the complete range of spectral types, including flat-spectrum sources (α≥ −0.5), steep-spectrum sources (α< −0.5), and peaked-spectrum sources. We show that there is no clear correlation between spectral type and size, but there is a correlation between the high-frequency spectral index and both object type and size. We also show that, to avoid biasing the data and to understand the various classes of active galactic nuclei (AGN) involved, the complete range of spectral types should be included in studying the general phenomenology of compact jetted AGN, and that complete samples must be used, selected over a wide range of frequencies. We discuss examples that demonstrate these points. We find that the high-frequency spectral indices of CSO-2s span −1.3 <α_hi< −0.3 and hence that radio spectral signatures cannot be used to discriminate definitively between CSO-2s, binary galactic nuclei, and millilensed objects, unless they haveα_hi> −0.3. 

The Carbon monOxide Mapping Array Project (COMAP) Pathfinder survey continues to demonstrate the feasibility of line-intensity mapping using high-redshift carbon monoxide (CO) line emission traced at cosmological scales. The latest COMAP Pathfinder power spectrum analysis is based on observations through the end of Season 2, covering the first three years of Pathfinder operations. We use our latest constraints on the CO(1–0) line-intensity power spectrum atz~ 3 to update corresponding constraints on the cosmological clustering of CO line emission and thus the cosmic molecular gas content at a key epoch of galaxy assembly. We first mirror the COMAP Early Science interpretation, considering how Season 2 results translate to limits on the shot noise power of CO fluctuations and the bias of CO emission as a tracer of the underlying dark matter distribution. The COMAP Season 2 results place the most stringent limits on the CO tracer bias to date, at ⟨T b⟩ < 4.8 μK, which translates to a molecular gas density upper limit ofρ_H2< 1.6 × 10⁸M_⊙Mpc⁻³atz~ 3 given additional model assumptions. These limits narrow the model space significantly compared to previous CO line-intensity mapping results while maintaining consistency with small-volume interferometric surveys of resolved line candidates. The results also express a weak preference for CO emission models used to guide fiducial forecasts from COMAP Early Science, including our data-driven priors. We also consider directly constraining a model of the halo–CO connection, and show qualitative hints of capturing the total contribution of faint CO emitters through the improved sensitivity of COMAP data. With continued observations and matching improvements in analysis, the COMAP Pathfinder remains on track for a detection of cosmological clustering of CO emission. 

The CO Mapping Array Project (COMAP) Pathfinder is performing line intensity mapping of CO emission to trace the distribution of unresolved galaxies at redshiftz ∼ 3. We present an improved version of the COMAP data processing pipeline and apply it to the first two Seasons of observations. This analysis improves on the COMAP Early Science (ES) results in several key aspects. On the observational side, all second season scans were made in constant-elevation mode, after noting that the previous Lissajous scans were associated with increased systematic errors; those scans accounted for 50% of the total Season 1 data volume. In addition, all new observations were restricted to an elevation range of 35–65 degrees to minimize sidelobe ground pickup. On the data processing side, more effective data cleaning in both the time and map domain allowed us to eliminate all data-driven power spectrum-based cuts. This increases the overall data retention and reduces the risk of signal subtraction bias. However, due to the increased sensitivity, two new pointing-correlated systematic errors have emerged, and we introduced a new map-domain PCA filter to suppress these errors. Subtracting only five out of 256 PCA modes, we find that the standard deviation of the cleaned maps decreases by 67% on large angular scales, and after applying this filter, the maps appear consistent with instrumental noise. Combining all of these improvements, we find that each hour of raw Season 2 observations yields on average 3.2 times more cleaned data compared to the ES analysis. Combining this with the increase in raw observational hours, the effective amount of data available for high-level analysis is a factor of eight higher than in the ES analysis. The resulting maps have reached an uncertainty of 25–50 μK per voxel, providing by far the strongest constraints on cosmological CO line emission published to date. 

We present updated constraints on the cosmological 3D power spectrum of carbon monoxide CO(1–0) emission in the redshift range 2.4–3.4. The constraints are derived from the two first seasons of Carbon monOxide Mapping Array Project (COMAP) Pathfinder line intensity mapping observations aiming to trace star formation during the epoch of galaxy assembly. These results improve on the previous Early Science results through both increased data volume and an improved data processing methodology. On the methodological side, we now perform cross-correlations between groups of detectors (“feed groups”), as opposed to cross-correlations between single feeds, and this new feed group pseudo power spectrum (FGPXS) is constructed to be more robust against systematic effects. In terms of data volume, the effective mapping speed is significantly increased due to an improved observational strategy as well as a better data selection methodology. The updated spherically and field-averaged FGPXS,C^~(k), is consistent with zero, at a probability-to-exceed of around 34%, with an excess of 2.7σin the most sensitive bin. Our power spectrum estimate is about an order of magnitude more sensitive in our six deepest bins across 0.09 Mpc⁻¹<k< 0.73 Mpc⁻¹, compared to the feed-feed pseudo power spectrum (FPXS) of COMAP ES. Each of these bins individually constrains the CO power spectrum tok P_CO(k) < 2400–4900 μK²Mpc²at 95% confidence. To monitor potential contamination from residual systematic effects, we analyzed a set of 312 difference-map null tests and found that these are consistent with the instrumental noise prediction. In sum, these results provide the strongest direct constraints on the cosmological 3D CO(1–0) power spectrum published to date. 

We present the first degree-scale tomography map of the dusty magnetized interstellar medium (ISM) from stellar polarimetry and distance measurements. We used the RoboPol polarimeter at Skinakas Observatory to conduct a survey of the polarization of starlight in a region of the sky of about four square degrees. We propose a Bayesian method to decompose the stellar-polarization source field along the distance to invert the three-dimensional (3D) volume occupied by the observed stars. We used this method to obtain the first 3D map of the dusty magnetized ISM. Specifically, we produced a tomography map of the orientation of the plane-of-sky component of the magnetic field threading the diffuse, dusty regions responsible for the stellar polarization. For the targeted region centered on Galactic coordinates (l,b) ≈ (103.3°, 22.3°), we identified several ISM clouds. Most of the lines of sight intersect more than one cloud. A very nearby component was detected in the foreground of a dominant component from which most of the polarization signal comes and which we identified as being an intersection of the wall of the Local Bubble and the Cepheus Flare. Farther clouds, with a distance of up to 2 kpc, were similarly detected. Some of them likely correspond to intermediate-velocity clouds seen in HIspectra in this region of the sky. We found that the orientation of the plane-of-sky component of the magnetic field changes along distance for most of the lines of sight. Our study demonstrates that starlight polarization data coupled to distance measures have the power to reveal the great complexity of the dusty magnetized ISM in 3D and, in particular, to provide local measurements of the plane-of-sky component of the magnetic field in dusty regions. This demonstrates that the inversion of large data volumes, as expected from the PASIPHAEsurvey, will provide the necessary means to move forward in the modeling of the Galactic magnetic field and of the dusty magnetized ISM as a contaminant in observations of the cosmic microwave background polarization. 

We present the first Bayesian method for tomographic decomposition of the plane-of-sky orientation of the magnetic field with the use of stellar polarimetry and distance. This standalone tomographic inversion method presents an important step forward in reconstructing the magnetized interstellar medium (ISM) in three dimensions within dusty regions. We develop a model in which the polarization signal from the magnetized and dusty ISM is described by thin layers at various distances, a working assumption which should be satisfied in small-angular circular apertures. Our modeling makes it possible to infer the mean polarization (amplitude and orientation) induced by individual dusty clouds and to account for the turbulence-induced scatter in a generic way. We present a likelihood function that explicitly accounts for uncertainties in polarization and parallax. We develop a framework for reconstructing the magnetized ISM through the maximization of the log-likelihood using a nested sampling method. We test our Bayesian inversion method on mock data, representative of the high Galactic latitude sky, taking into account realistic uncertainties from Gaia and as expected for the optical polarization survey P ASIPHAE according to the currently planned observing strategy. We demonstrate that our method is effective at recovering the cloud properties as soon as the polarization induced by a cloud to its background stars is higher than ~0.1% for the adopted survey exposure time and level of systematic uncertainty. The larger the induced polarization is, the better the method’s performance, and the lower the number of required stars. Our method makes it possible to recover not only the mean polarization properties but also to characterize the intrinsic scatter, thus creating new ways to characterize ISM turbulence and the magnetic field strength. Finally, we apply our method to an existing data set of starlight polarization with known line-of-sight decomposition, demonstrating agreement with previous results and an improved quantification of uncertainties in cloud properties. 

Context. Optical polarimeters are typically calibrated using measurements of stars with known and stable polarization parameters. However, there is a lack of such stars available across the sky. Many of the currently available standards are not suitable for medium and large telescopes due to their high brightness. Moreover, as we find, some of the polarimetric standards used are in fact variable or have polarization parameters that differ from their cataloged values. Aims. Our goal is to establish a sample of stable standards suitable for calibrating linear optical polarimeters with an accuracy down to 10 −3 in fractional polarization. Methods. For 4 yr, we have been running a monitoring campaign of a sample of standard candidates comprised of 107 stars distributed across the northern sky. We analyzed the variability of the linear polarization of these stars, taking into account the non-Gaussian nature of fractional polarization measurements. For a subsample of nine stars, we also performed multiband polarization measurements. Results. We created a new catalog of 65 stars (see Table 2) that are stable, have small uncertainties of measured polarimetric parameters, and can be used as calibrators of polarimeters at medium and large telescopes. 

Context.Blazars are beamed active galactic nuclei (AGNs) known for their strong multi-wavelength variability on timescales ranging from years down to minutes. Many different models have been proposed to explain this variability. Aims.We aim to investigate the suitability of the twisting jet model presented in previous works to explain the multi-wavelength behaviour of BL Lacertae, the prototype of one of the blazar classes. According to this model, the jet is inhomogeneous, curved, and twisting, and the long-term variability is due to changes in the Doppler factor due to variations in the orientation of the jet-emitting regions. Methods.We analysed optical data of the source obtained during monitoring campaigns organised by the Whole Earth Blazar Telescope (WEBT) in 2019–2022, together with radio data from the WEBT and other teams, andγ-ray data from theFermisatellite. In this period, BL Lacertae underwent an extraordinary activity phase, reaching its historical optical andγ-ray brightness maxima. Results.The application of the twisting jet model to the source light curves allows us to infer the wiggling motion of the optical, radio, andγ-ray jet-emitting regions. The optical-radio correlation shows that the changes in the radio viewing angle follow those in the optical viewing angle by about 120 days, and it suggests that the jet is composed of plasma filaments, which is in agreement with some radio high-resolution observations of other sources. Theγ-ray emitting region is found to be co-spatial with the optical one, and the analysis of theγ-optical correlation is consistent with both the geometric interpretation and a synchrotron self-Compton (SSC) origin of the high-energy photons. Conclusions.We propose a geometric scenario where the jet is made up of a pair of emitting plasma filaments in a sort of double-helix curved rotating structure, whose wiggling motion produces changes in the Doppler beaming and can thus explain the observed multi-wavelength long-term variability.